Shell model studies of the proton drip line nucleus 106 Sb 
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We present results of shell model calculations for the proton drip line nucleus 106 Sb. The shell 
model calculations were performed based on an effective interaction for the 2sld0g 7 / 2 0hlln/2 shells 
employing modern models for the nucleon-nucleon interaction. The results are compared with the 
recently proposed experimental yrast states. A good agreement with experiment is found lending 
support to the experimental spin assignements. 
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Considerable attention is at present being devoted to 
the experimental and theoretical study of nuclei close to 
the limits of stability. Recently, heavy neutron deficient 
nuclei in the mass regions of A = 100 have been studied, 
and nuclei like 100 Sn and neighboring isotopes have been 
identified Moreover, the proton drip line has been 

established in the A = 100 and A = 150 regions PJ and 
nuclei like 105 Sb and 109 I have recently been established 
as ground-state proton emitters (^||. The next to drip 
line nucleus for the antimony isotopes, 106 Sb with a pro- 
ton separation energy of ~ 400 keV, was studied recently 
in two experiments and a level scheme for the yrast states 
was proposed in Ref. 0. 

The aim of this work is thus to see whether shell-model 
calculations, which employ realistic effective interactions 
based on state of the art models for the nucleon-nucleon 
interaction, are capable of reproducing the experimental 
results for systems close to the stability line. Before we 
present our results, we will briefly review our theoretical 
framework. In addition, we present results for effective 
proton and neutron charges based on perturbative many- 
body methods. These effective charges will in turn be 
used in a shell-model analysis of E2 transitions. 

The aim of microscopic nuclear structure calculations 
is to derive various properties of finite nuclei from the un- 
derlying hamiltonian describing the interaction between 
nucleons. We derive an appropriate effective two-body 
interaction for valence neutrons and protons in the single- 
particle orbits 2si/2, W5/2, W3/2, O57/2 and Ohn/2- As 
closed shell core we use 100 Sn. This effective two-particle 
interaction is in turn used in the shell model model study 
of 106 Sb. The shell model problem requires the solution 
of a real symmetric n x n matrix eigenvalue equation 



H |tf fc ) = E k |tf fc ) , 
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with k — 1, . . . , K. At present our basic approach to 
finding solutions to Eq. (|l|) is the Lanczos algorithm, an 
iterative method which gives the solution of the lowest 
eigenstates. The technique is described in detail in Ref. 
|f, see also Ref. §. 

To derive the effective interaction, we employ a pertur- 
bative many-body scheme starting with the free nucleon- 
nucleon interaction. This interaction is in turn renor- 
malized taking into account the specific nuclear medium. 



The medium rcnormalized potential, the so-called G- 
matrix, is then employed in a perturbative many-body 
scheme, as detailed in Ref. (h]] and reviewed briefly be- 
low. The bare nucleon-nucleon interaction we use is the 
charge-dependent meson-exchange model of Machleidt 
and co-workers jll), the so-called CD-Bonn model. The 
potential model of Ref. jllj is an extension of the one- 
boson-exchange models of the Bonn group fl2|| , where 
mesons like 7r, p, 77, 5, uj and the fictitious a meson are 
included. In the charge-dependent version of Ref. |fll| , 
the first five mesons have the same set of parameters for 
all partial waves, whereas the parameters of the a meson 
are allowed to vary. 

The first step in our perturbative many-body scheme is 
to handle the fact that the repulsive core of the nucleon- 
nucleon potential V is unsuitable for perturbative ap- 
proaches. This problem is overcome by introducing the 
reaction matrix G given by the solution of the Bcthc- 
Goldstone equation 



G = V + V- 
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where uj is the unperturbed energy of the interacting nu- 
cleons, and Ho is the unperturbed hamiltonian. The op- 
erator Q, commonly referred to as the Pauli operator, is 
a projection operator which prevents the interacting nu- 
cleons from scattering into states occupied by other nu- 
cleons. In this work we solve the Bethe-Goldstone equa- 
tion for five starting energies to, by way of the so-called 
double-partitioning scheme discussed in e.g., Ref. [ fTpf . 
A harmonic-oscillator basis was chosen for the single- 
particle wave functions, with an oscillator energy TiQ, 
given by htt = 45^" 1 / 3 - 25A- 2 / 3 = 8.5 MeV, A = 100 
being the mass number. 

Finally, we briefly sketch how to calculate an effec- 
tive two-body interaction for the chosen model space in 
terms of the G-matrix. Since the G-matrix represents 
just the summmation to all orders of ladder diagrams 
with particle-particle intermediate states, there are obvi- 
ously other terms which need to be included in an effec- 
tive interaction. Long-range effects represented by core- 
polarization terms are also needed. The first step then is 
to define the so-called Q-box given by 
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The Q-box is made up of non-folded diagrams which are 
irreducible and valence linked. The projection operators 
P and Q define the model space and the excluded space, 
respectively, with P + Q = I. All non-folded diagrams 
through third order in the interaction G are included in 
the definition of the Q-box while so-called folded dia- 
grams are included to infinite order through the summa- 
tion scheme discussed in Refs. Jl0| , |l3[ ] . 

Effective interactions based on the CD-Bonn nucleon- 
nucleon interaction have been used by us for several mass 
regions, and give in general a very good agreement with 
the data, see Refs. Jlf^l. 

In addition to deriving an effective interaction for the 
shell model, we present also effective proton and neutron 
charges based on our perturbative many-body methods. 
These charges are used in our studies of available E2 data 
below. In this way, degrees of freedom not accounted for 
by the shell-model space are partly included through the 
introduction of an effective charge. The effective single- 
particle operators for the effective charge are calculated 
along the same lines as the effective interaction. In nu- 
clear transitions, the quantity of interest is the transition 
matrix element between an initial state and a final 
state 1 4"/) of an operator O (here it is the E2 operator) 
defined as 



o fi = 



(4) 



Since we perform our calculation in a reduced space, the 
exact wave functions \& 1 j are not known, only their pro- 
jections onto the model space. We are then con- 
fronted with the problem of evaluating Ofi when only 
the model space wave functions are known. In treating 
this problem, it is usual to introduce an effective opera- 
tor O e g different from the original operator O defined by 
requiring 



O fi = <*/|Oeff|$i> 



(5) 



The standard scheme is then to employ a perturbative ex- 
pansion for the effective operator, see e.g. Refs. (l^ |2(| . 

To obtain effective charges, we evaluate all effec- 
tive operator diagrams through second-order, excluding 
Hartree-Fock insertions, in the G-matrix obtained with 
the CD-Bonn interaction. Such diagrams are discussed 
in the reviews by Towner |18) and Ellis and Osnes [fl9|| . 
The state dependent effective charges are listed in Table 
H for the diagonal contributions only. In order to repro- 
duce the experimental B(E2; 4^~ —^2^) transition of Ref. 
0, the authors introduced effective charges e p — 1.72e 
and e n = 1.44e for protons and neutrons, respectively. 
We see from Table | that the microscopically calculated 
values differ significantly from the above values from Ref. 
M. This could however very well be an artefact of the 



chosen model space and effective interaction employed in 
the shell model analysis of Ref. 0. The reader should 
also keep in mind that our model for the single-particle 
wave functions, namely the harmonic oscillator, may not 
be the most appropriate for the proton single-particle 
states, since the proton separation energy is of the or- 
der of some few keV. When compared with the theoret- 
ical calculation of Sagawa et al. f2l|| , our neutron effec- 
tive charges agree well with theirs, whereas the proton 
effective charge deduced in Ref. [£fj is slightly larger, 
e p ~ 1.4e. We note also that in the Hartree-Fock cal- 
culation with a Skyrme interaction and accounting for 
effects from the continuum, Hamamoto and Sagawa p^j 
obtained effective charges of e„ = 1.35e and e p — l.Oe for 
100 Sn. Below we will allow the effective charges to vary 
in order to reproduce as far as possible the experimental 
value of 2.8(3) W.u. for the transition B{E2;A J [ -> 2^). 
There we will also relate the theoretical values for the 
effective charges to those extracted from data around 
A = 100, see e.g., Ref. |§. 

The calculations were performed with two possible 
model spaces, one which comprises all single-particle or- 
bitals of the lrf5/20g7/2lc?3/22si/20/2.ii/2 shell and one 
which excludes the 0hn/ 2 orbit. The latter model space 
was employed by the authors of Ref. |?J in their shell 
model studies. Since the single-neutron and single-proton 
energies with respect to 100 Sn are not well-established, 
we have adopted for neutrons the same single-particle 
energies as used in large-scale shell-model calculations 
of the Sn isotopes, see Refs. The neutron single- 
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= 2.45 MeV and 



£0/111/2 — e ic*5/2 = 3-2 MeV. These energies, when em- 
ployed with our effective interaction described above, 
gave excellent results for both even and odd tin isotopes 
from 102 Sn to 116 Sn. The proton single-particle ener- 
gies are less established and we simply adopt those for 
the neutrons. Since the proton separation energy is of 
the order of ~ 400 kev, it should suffice to carry out a 
shell-model calculation with just the 1^5/2037/2 orbits for 
protons. The total wave function, see the discussion be- 
low, is however to a large extent dominated by the 1^5/2 
orbital for protons, with small admixtures from the 0^7/2 
proton orbital. The influence from the other proton or- 
bits is thus minimal. 

The resulting eigenvalues are displayed in Fig. [I] for 
the two choices of model space together with the ex- 
perimental levels reported in Ref. 0. Not all exper- 
imental levels have been given a spin assignment and 
all experimental spin values are tentative. The label 
FULL stands for the model space which includes all or- 
bits from the lrf5/20g7/2ld 3 /22s 1 / 2 0/i 11 /2 shell while RE- 
DUCED stands for the model space where the O/in/2 
orbit has been omitted. As can be seen from Fig. 
[I], the agreement with experiment is also rather good, 
with the model space which includes all orbitals of the 
1^5/2057/2lrf3/22si/20/i n / 2 shells being closest to the ex- 
perimental level assignements. The reader should how- 
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ever note that in Ref. j7| it is not excluded that the 
ground state could have spin 1 + , which means that the 
experimental spin values in Fig. |l| should be reduced by 
1. In our theoretical calculations we obtain in addition 
to a state with spin 1 + , also a state with 3 + not seen 
in the experiment of Ref. 0]. In case the ground state 
turns out to have spin 1 + , the reduced model space in 
our calculations would yield a better agreement with the 
data. 

The wave functions for the various states are to a large 
extent dominated by the 0<? 7 / 2 and ldg/2 single-particle 
orbits for neutrons (y) and the lds/2 single-particle orbit 
for protons (ir). The h'Ogi/2 and ^1^/2 single-particle 
orbits represent in general more than ~ 90% of the to- 
tal neutron single-particle occupancy, while the 
single-particle orbits stands for ~ 80 — 90% of the pro- 
ton single-particle occupancy. The other single-particle 
orbits play an almost negligible role in the structure 
of the wave functions. The only notable exception is 
the 7 1 state where 7r0<77/2 stands for the 84% of the 
proton single-particle occupancy. This has also impor- 
tant repercussions on the contributions to the measured 
E2 transition £>(_E2;4^ — » 2f), where the stucture of 
the wave functions of the 4^ and 2^ states are to a 
large extent dominated by the 1^1^5/2 and 7rld5/2 single- 
particle orbits. The O57/2 orbits play a less significant 
role in the structure of the wave functions, and since 
the Ogj/2 W5/2 transition matrix element tends to be 



weaker than the one between left 
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sition will be dominated by the latter contributions. As 
also noted by Sohler et al. |7j], the E2 transition is dom- 
inated by neutron contributions. This can also be seen 
from Fig. || where we show the result for the above E2 
transition as function of different choices for the effective 
charges. We see that the largest change in the value of 
the E2 transition takes place when we vary the effective 
charge of the neutron, whereas when the proton charge 
is changed, the percentual change is smaller. Further- 
more, if we use the largest values for effective charges of 
Table H, namely e n = 0.72e and e p = 1.16e, we obtain 
1.84 W.u. for the E2 transition. Compared with the ex- 
perimental value of 2.8(3) W.u. this may indicate that 
both the proton and the neutron effective charges should 
be slightly increased. From Fig. || we see that effective 
charges of e„ = 0.9e and e p = 1.4e ± 0.2e seem to yield 
the best agreement with experiment, although neutron 
charges of e n = 0.8e and e n = l.Oe yield results within 
the experimental window of Fig. |^ . The neutron ef- 
fective charges would agree partly with those exctracted 
from the data in the Sn isotopes 21 and from theoret- 
ical calculations of E2 transitions in heavy Sn isotopes 
fusf , where a value e n ~ 1 is adopted in order to repro- 
duce the data. The calculated effective charges of Table 
H are however on the lower side. However, the deduced 
effective charges from the B(E2; 6^ — > 4^) transitions in 
102 Sn H and 104 Sn @ indicate that e„ ~ 1.6 - 2.3e, 
depending on the effective interaction employed in the 
shell-model analyses. Clearly, the effective interaction 



which is used, and its pertinent model space, approx- 
imations made in the many-body formalism etc., will 
influence the extraction of effective charges. This no- 
table difference in the effective charges could be due to 
the fact that the B(E2\6f -> 4f) transitions in 102 Sn 
and 104 Sn involve configurations not accounted for by 
the 1^5/20(77/21^3/22,31/20/111/2 model space. A proton 
effective charge of e p ~ 1.4e is close to values inferred 
from experiment for the N — 50 isotones, see Ref. jjJU, 
shell-model calculations of E2 transitions for the N = 82 
isotones 25 1 and the theoretical estimates of Ref. |2lj] . 

In summary, a shell-model calculation with realistic ef- 
fective interactions of the newly reported low-lying yrast 
states of the proton drip line nucleus 106 Sb, reproduces 
well the experimental data. Since the wave functions of 
the various states are to a large extent dominated by neu- 
tronic degrees of freedom and neutrons are well bound 
with a separation energy of ~ 8 MeV, this may explain 
why a shell-model calculation, within a restricted model 
space for a system close to the proton drip line, gives 
a satisfactory agreement with the data. In order to re- 
produce the experimental B(E2;4^ — > 2f) transition, 
we obtained effective charges from our shell-model wave 
functions of e n = 0.8 — l.Oe and e p = 1.4e ± 0.2e. Our 
microscopically calculated effective charges are however 
slightly smaller, e n — 0.5 — 0.7e and e p = 1.1 — 1.2e. 
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TABLE I. Proton and neutron effective charges relative to 



Sn for the ld 5/2 , 0g r/2 , ld 3 / 2 , 2s 1/2 and 0h 11/2 single par- 
ticle orbitals. 





Proton 


Neutron 


1^5/2 


1.06e 


0.53e 


0^7/2 


1.15e 


0.72e 


W3/2 


1.04e 


0.52e 


0^11/2 


1.16e 


0.51e 



4 



MeV 



4 - 



(7+) 



13+ 



12+ 



11+ 



13+ 



12+ 



10+ 11+ 

9+ 

(10+) 10+ 



9+ 

8+ 

(8+) 8+ 

7+ 



7+ 



6+ 



6+ 

5+ 5+ 

3+ 3+ 

4+ 4+ 

4+) 1+ 2+ 

2+) 2+ 1+ 

EXPT FULL REDUCED 



FIG. 1. Low-lying states for Sb, theory and experiment. Energies in MeV. FULL means the model space which comprises 
all single-particle orbits, 2sld0g 7 / 2 0hllii/ 2 - REDUCED means that the 0/illn/ 2 single-particle orbit is not included. 
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FIG. 2. Value for the B{E2\A J [ —* 2~l) transition as function of different effective charges in units of W.u. The horizontal 
lines represent the experimental window, with a value of 2.8 ± 0.3 W.u. 
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